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Abstract 

The character of the interaction between an impurity vibrational mode and a heat 
bath leads to certain peculiarities in the relaxational dynamics of the excited states. 
We derive a non-Markovian equation of motion for the reduced density matrix of 
this system which is valid for initial, intermediate and kinetic stages of the relax- 
ation. The linear phase-sensitive character of the interaction ensures the ultrafast 
disappearance of the quantum interference of the initially superpositional state and 
the effect of classical squeezing of the initially coherent state. On the other hand, the 
second power interaction induces a partial conservation of the quantum interference. 
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Motivation and model 

Time-resolved experimental techniques [1] sometimes allow to detect quantum su- 
perpositional states in different physical systems [2,3]. Such states provide a basis 
for the possibility of applications like quantum computers [4]. Unfortunately, these 
interesting states inevitably decay due to the coupling with a heat bath containing 
many degrees of freedom [5] . Below we consider how the character of the coupling 
with the bath influences the dynamics of the superpositional state of the impurity 
vibrational mode. Describing this mode as a harmonic oscillator, one writes the 
interaction part of the Hamiltonian as 

H: = hY,K{uj^) {b+f{a,a+) + bj+{a,a+)) , (1) 

where the function K{u!^) describes the intensity of the interaction between the 
vibronic mode and the bath mode operators b^. f{a, a+) is a function of the vibronic 
mode operators. We consider two cases, namely f{a, a"*") — a + a"^, corresponding 
to the case of the linear phase-sensitive interaction, and f{a,a~^) = yielding 
a quadratic interaction in the rotating wave approximation. Below we are using 
"single" and "double" to indicate the baths with these two types of interactions. 

Single bath 

The first case describes the processes when system and bath are exchanging 
one quantum. Such a behavior is provided by the physical situation when the 
majority of the bath modes contributing to dephasing and thermalization of the 
system have approximately the same frequency lo^ as the system mode uj. Applying 
the formalism of the time evolution operator, and restricting to the second order 
cumulant expansion [6], we obtain the non-Markovian master equation [7] for the 
reduced density matrix o of the impurity vibrational mode: 
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— = —iLj[a'^a, a 



(2) 



+ (7„+i + 7*) [aa, a+ + a] 



+ (7^+1 +7n)[a+ + a,aa+] 



+ (7n+i + 7n)[«^cr'«^ + «] 



+ (7*+i +7n)[a+ + a,cTa]. 



We obtain four relaxation functions 7, which are shown in Fig. la. They originate 
from the correlations between the operators a{t) and a~^{t) of the system mode uu 
and the memory kernels of the bath, like ^6^(0)6^/(t)^, which appear in the second 
order cumulant expansion of the evolution operator. These coefficients are found to 
be time-dependent, jn+i and 7„ describe the situation when emission-like processes 
7„+i prevail over the absorption-like processes 7„, where n denotes the number of 
quanta in the bath modes. The functions 7n+i, 7n correspond to the reverse situation 
and are always small. 

The evolution of the different initial states of the system was found [6] in ana- 
lytical form for two cases, namely: the initial stage of relaxation, when all of the 
relaxation functions are linear in time, and for the kinetic stage of the relaxation, 
when the coefficients %, 7„+i vanish, while 7^, 7„+i become constants Fi^i and 
ri(ni -|- 1), respectively. Here rii — n{cu) indicates the number of quanta in the 
bath mode at the system frequency. The analytical solution is based on the gener- 
ating function formalism [6] . The approaches corresponding to the above stages are 
applied to the evolution of the superposition of two coherent states |a) and e"^ |— «) 



where is the normalization constant, \a) is obtained by displacement of the vac- 
uum state |0) by a as \a) = exp(aa+ — a*a) |0). The dependence of the probability 




(3) 
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density P on coordinate Q and time is found to consist of classical and interference 
parts [7]: 

P{Q, t) = Pclass(Q, t) + PintiQ, t). (4) 

The interference part behaves somewhat different under each approach, compare Fig. 
lb. Nevertheless, energy relaxation of Pciass and decoherence of Pint occur on different 
time scales, independent on the approach applied: The quantum interference Pint 
disappears faster. The phase-sensitive character of the relaxation induces small 
oscillations of the broadening of the initially coherent wave packet [7] . 



Double bath 

When the system interacts with a bath having the maximum of its mode density 
at twice the oscillator frequency, then processes occur in which the system loses 2 
quanta and the bath obtains one quantum. The reverse processes are also allowed. 
To describe such a behavior we use /(a, a"*") = in the interaction Hamiltonian 
(1). The kinetic stage of the evolution of such a system follows the master equation 



'dt 



-tUJ 



a a, a 
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+ r2(n2 + l){ 
+ r2n2{[(a+)V, 



aV, (a+)2 



+ 
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where r2 = T:K'^g2 is the decay rate of the vibrational amplitude. Here, the number 

of quanta in the bath mode n2 = n{2uj), the coupling function K = K{2lj), and 
the density of bath states g2 — g{2uj) are evaluated at the double frequency of the 
selected oscillator. 

Rewritten in the basis of eigenstates \n) of the unperturbed oscillator this mas- 
ter equation contains only linear combinations of such terms as (Tm,n — (m| cr 



Cm+2,n+2, and CFjn-2,n-2- It ensures, in effect, a possibility to distinguisfi even and 
odd initial states of the system. The odd excited state |1) cannot relax into the 
ground state |0), but the even excited state |2) can. 

The evolution of the system after preparation under different initial conditions 
was simulated numerically. The equations of motion of the density matrix elements 
are integrated using a fourth order Runge-Kutta algorithm with stepsize control. To 
make the set of differential equations a finite one we restrict the number of levels by 
m,n < 20. 

One of the representative examples is the behavior of the initially coherent state 
of the system. For comparison with the usual behavior we have used the time 
dependence of the mean value of the coordinate. One can see in Fig. 2a that the 
relaxation consists of two stages. The coordinate mean value of the usual system 
decreases with a constant rate during both stages. The same initial value of the 
system coupled to a double bath shows a fast decrement in the first stage and 
almost no decrement in the second. 

The question about the influence of the type of the bath on the evolution of 
the superpositional state is of special importance. The simulation was made at 
a temperature of UbT = 2/ia;/ln3, corresponding to n(2a;) = 0.5. To allow the 
comparison, we have simulated the evolution of the same superpositional states 
twice. The only difference was the type of the baths, namely single and double 
bath. The coordinate representation of the wave packets is presented in Fig. 2b 
and 2c. The same value of the relaxation Fi = F2 provides, however, different 
results. In the system coupled to the single bath the quantum interference disappears 
already during the first period, while the amplitude decreases only slightly. In the 
opposite case, the system coupled to the double bath leaves the interference almost 
unchanged, although a fast reduction of the amplitude occurs. 
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Therefore, the second system partially conserves the quantum superpositional 
state. The experimental investigation of systems displaying such properties [2] is 
necessary, both for extracting typical parameter ranges for theoretical models and 
for practical applications like quantum computation and quantum cryptography. 

Conclusions 

The character of the coupling with a heat bath plays a dominating role in the time 
evolution of different excited states of a vibrational mode. Linear coupling ensures 
the ultrafast decay of the interference part of the superpositional state. This result 
remains true even if different approaches are applied. The quadratic type of coupling 
gives the same time scales both for the amplitude and interference relaxations. Some 
initial coherent properties, like the distinction between even and odd levels, survive 
for long time scales. 
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FIGURES 

FIG. 1. (a) Time-dependent relaxation coefficients jn+i (soUd line), jn (diamonds), 
7„+i (boxes), 7„ (crosses) of the master Eq. (2), calculated for a bath containing 60 modes 
in the range [0,6a;], with coupling function K{u>^) = 1. (b) Time dependence of Pint 
from Eq. (4), when the heat bath can be described with only two constants Fi = 0.25a;, 
ni = 0.4 for a; = 1, Q{t = 0) = A, (p = 0. Crosses: derivation valid for the initial 
stage of the relaxation, solid line: result for the kinetic stage of the relaxation, diamonds: 
corresponding result obtained for the linear interaction in the rotating wave approximation 
/(a,a+) = a. 

FIG. 2. Difference between double and single baths, (a) Evolution of the coherent 
state: Mean value of the coordinate, for Q{t = 0) = —2.2. Diamonds: system coupled to 
the single bath, Fi = 0.15a;. Solid line: system coupled to the double bath, F2 = 0.5a;. 
(b) Time evolution of the superposition of coherent states with initial separation 2Qq = 8 
interacting with single bath (b) for a; = 1, Fi = 0.005a;, n{uj) = 1.36. (c) Same as (b), but 
for a double bath with F2 = 0.005a;, n(2a;) = 0.5. 
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